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ABSTRACT 



We calculate the expected number of type la supernovae (SN la) in the core- 
degenerate (CD) scenario and find it to match observations within the uncertainties 
of the code. In the CD scenario the super-Chandrasekhar mass white dwarf (WD) 
is formed at the termination of the common envelope phase from a merger of a WD 
companion with the hot core of a massive asymptotic giant branch (AGB) star. We 
use a simple population synthesis code that avoids the large uncertainties involved in 
estimating the final orbital separation of the common envelope evolution. Instead, we 
assume that systems where the core of the secondary AGB star is more massive than 
the WD remnant of the primary star merge at the termination of the common envelope 
phase. We also use a simple prescription to count systems that have strong interaction 
during the AGB phase, but not during the earlier red giant branch (RGB) phase. That 
a very simple population synthesis code that uses the basics of stellar evolution ingre- 
dients can match the observed rate of SN la might suggest that the CD-scenario plays 
a major role in forming SN la. 



INTRODUCTION 



There are four basic theoretical scenarios for the formation of the progenitors of Type la 
supernovae (SNe la). The goal of these scenarios is to form a carbon-oxygen (CO) white dwarf 
(WD) with a mass close to and ab ove the Chandrasekha r mass limit Mch- Such WDs will go 
through a thermonuclear detonation (Hovle & Fowler I96 0]) that is observed as an SN la. (i) In the 
singl e degenerate (SD) scenario (e.g., Whelan &z Ibenl ' l973 : Nomoto 1982 : Han &: Podsiadlowski 



2004) a WD grows in mass through accretion from a no n-degenerate stella r companion. However, 
it seems that the WD-mass increase is very limited, e.g., Ildan et al.l (|2012l ) for a recent paper and 
references therein. Ruiter et al. (2011) consider the helium-rich donor scenario (HeRS) to be in 
a category separate from the canonical SD scenario. In the HeRS a He-rich star, degenerate or 
not degenerate, tran sfers mass to a CO white dwarf which explodes as it approaches the Chan- 
drasekhar mass (e.g. Ilben et al.lll987l ). Like the SD scenario, the delay time is strongly dependent 
on evolutionary timescales of the stars. (ii) In the double degenerate (DD) scenario (jWebbink 
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1984 ; llben fc Tutukovlll984l ; see Ivan Kerkwiik et al.l I201Q for a paper on sub-Chandrasekhar mass 



remnants) two WDs merge after losing energy a nd angular momentum through the radiation of 



gravitational waves (jTutukov &: Yungelsonlll979l ). (Hi) In the core-degenerate (CD) scenario for 



the formation of SN la the Chandrasekhar or super-Chandrasekhar mass WD is formed at the 
termination of the common envelope (CE) phase or during the planetary nebula phase, from a 



( Kashi Soker 



2011 



Ilkov &; Sokei 


2012; 


Soker 


2011 



Soker et al. 



2012). Observations and theo- 



retical studies cannot teach us yet w hether all scenario, or only one or two of these can work (e.g., 



Livio 



2001 



Maoz 



2010 



Howeilll201ll ), (iv) It is likely th at a sub-Chandrasekhar mass WD can ex 



plode as a SN la via the 'double-detonation' mechanism (jWooslev fc Weaver 



1994 



Livne Arnett 



19951 ). In this model, a sub-Chandrasekhar mass WD accumulates a layer of helium-rich material 
on the surface, which under the ri ght c onditions can detonate, leading to a detonation near the 



center of the CO WD (|Fink et al.ll2010l ). Ruiter et al. (2011) performed a population synthesis 
study of the double-detonation scenario considering helium-rich donors, and found a bimodal delay 
time distribution, the shape depending on whether non-degenerate donors (< 1 Gyr) or degenerate 
donors (> 1 Gyr) dominate. Regarding both the DD and the double-detonation scenarios, we note 
that sub-Chandrasekhar SN la explosions are s upported also by the population of WD binaries in 
the solar neighborhood (jBadenes &: Maoal2012l ). 



in the Dast ( Soarks & Stecher 


1974; 


Livio & Riessl ( 


2003|) suggested that 



the merger of the WD with the AGB core leads to a SN la that occurs at the end of the CE phase or 
sh ortly after, and can explain the presence of hydrogen lines. This idea was further studied recently 
bv lSoker et al.1 ( 20120. In the CD scenari o the possibility of a very long time delay (up to 10 10 yr) 



is considered as well (lllkov &: Sokerl l2012). However, if that delay mechanism does not wor k, then 
the CD scenario ca n explain SN la only in star forming galaxies. iMennekens et al.l (|2010i ) and 



Ruiter et al.l (|2008l ) find in their population synthesis calculations that many systems end in a core- 



Anandlll965l ; lOstriker &: Bodenheimerlll968l ; Uenishi et al 



Uenishi et al. 


2003 




Yoon & L anger 


2005), and pos- 


■ 


Kundu & Mukhopadhvav 


2012; 


Garcfa-Berro et al. 



2012 ). the super-Chandrasekhar WD does not explode. 



There are two key ingredients in the CD scenario, in addition to the common condition that 
the remnant mass be > 1.4Mq. (1) The merger should occur while the core is still l arge, hence hot. 
This l imits the merger to occur within ~ 10 5 yr after the common envelope phase. iKashi fc Soker 
(|201ll ) showed that this condition can be met whe n the AGB star is massive, and some of the 
ejected CE gas might fall back (see also |SokerJ|2012j). (2) The delay between merger and explosion 
should be up to ~ 10 10 yr if this scenario is to account for SNe la in old-stellar populations, in 
addition to SNe la in young stellar populations. WDs with delay times have masses in the range 
~ 1.4 — 1.5M , and indeed can have a long delay (Jorge Rueda, 2012, private communication). 
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In previous papers another ingredient was used. It asserted that the hot core is more massive than 
the compani on cold WD , This is true in most cases, but not in all. In some rare cases this not need 
be the case (jSoker et al.ll2012l ). and there is a short delay to explosion. 



These ingredients involve physical processes much different from those in the DD scenario, and 
hence make the CD a distinguished scenario, rather than a branch of the DP scenario. The merger 
of the core while it is still hot might prevent an early ignit ion of carbon (lYoon et al.l 120071 ). which 
is one of the theoretical problems of the DD scenario (e.g., ISaio &: Nomotoll2004l ). As for the delay 
time, in t he CD scenario it is due to the spinning-down time of the WD set by magneto-dipole 
radiation (lllkov <fe Sokerll2012l ). while in the DD scenario the delay is due to the spiraling-in time 
due to gravitational radiation. 



Stabilizing rapidly rotating super-Chand rase khar WDs is a delicate matter (e.g.. lBoshkavev et al. 



20121 . lYoon & Langerll20ol . Ichen k Lilbood and lHachisu et al.ll2012al lb1 in the SD degenerate 



sce- 
nario). The strong magnetic fields required in the present model most likely will enforce a rigid 
rotation within a short time scal e due the WD being a perfect conductor. The critical mass of 



rigidly rotating WDs is 1.48M0 (jYoon &: Langerl I2004J and references therein). This implies that 



WDs more massive than 1.48Mq will explode in a relatively short time. The similarity of most SN 
la suggests that their progenitors indeed come from a narrow mass range. This is ~ 1.4 — 1.48Mq 
in the CD scenario. This property of the magneto-dipole radiation torque spinning-down mecha- 
nism, that only WDs with Mwd ^ 1.48Mq can slow down on a very long time scale (e.g., Jorge 



less luminous (e.g.. iHowell 



Howell 


2001 


Smith et al. 


2011 



and that very massi ye SN la progenito rs occur 



in galaxies where star formation is expected, i.e., spirals and irregulars (jScafzo et al.ll2012l ). 



In the present paper we estimate the number of SNe la from the CD scenario. Because of 
several large uncertainties, we perform a crude estimate of the SN la progenitor formation rate. In 
section [2] we point to an interesting coincident of the Chandrasekhar mass with another mass in 
the CD and DD scenarios. The assumptions and set up of the population synthesis are given in 
section [3j The results on the formation rate are given section U together with a comparison with 
observations, including with core-collapse supernovae (CCSNe) rate. Our summary is in section [5j 



2. THE ~ 1.35M COINCIDENCE 

The electron degeneracy pressure operating in two different stellar phases leads to an interesting 
coincidence in the frame of the scenario discussed here. Relativistic electron degeneracy pressure 
sets the Chandrasekhar mass limit of Mch = 1.4M0. Non-relativistic electron degeneracy pressure 
operates in the core of low mass red giant branch (RGB) stars. These stars reach a relatively 
large radius on the RGB. The combined mass of the two WD remnants of a binary system where 
both stars develop a degenerate core on the RGB is limited by ~ 1.3M0. When evolutionary 
considerations are added (see below) this mass is very close to Mch- 
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Stars that have a main sequence (MS) mass of Mms i$ 2.3M0 develop a degenerate core 
and reach a relatively large radius on the RGB, not much smaller than the maximum radius they 
achieve on the asymptotic giant branch (AGB). Stars with Mms ^ 2.3M , on the other hand, reach 
a much larger radius on the AGB relative to the radius they achieve on the RG B. The following 
approximation (up to 20%) can be fitted to the results of llben h Tutukovl (|1985l ) (their fig. 31) for 
the ratio of the maximum radius on the AGB to that on the RGB as function of the MS mass 



\og(R A /R R ) = 3.71og 2 (M/M ) - O.371og(M/M ) + 0.16, M < 2.25M G 



and 



log(R A /R R ) = 2.2 - 1.81og(M/M e 



M > 2.35Mr 



(1) 



(2) 



The results of I lb en Tutukovl ()1985l ) together with the fitting formula are plotted in Figure [TJ 
This implies that for stars with Mms ^ 2.3M in binary systems there is a large range of orbital 
separations where strong binary interaction occurs during the AGB phase but not during the RGB 
phase of the star. This condition is required to form a CO WD merger remnant. 



It so happens that a star of MS mass of Mms 



2.3 that evolves witho ut any disturbances 

This leads to the 



Catalan et al 



2008 



2009). 



forms a WD of mass M W d(2.3) ~ 0.65M Q (e.g. 
following argument. Consider binary systems where for both stars there is a large range of the 
binary orbital separation for a strong binary interaction to occur on the AGB, but not on the 
RGB. Namely, both stars have Mms > 2.3M when they leave the MS. The secondary can start 
the MS with Mm so < 2.3M , but then accretes mass from the primary star to the point where 
its mass is M2 new > 2.3M . We use the condition that after mass transfer the secondary mass is 
M2 n ew > 2.3M . We also use the condition M2 > 1M on the initial mass of the secondary star. 
This ensures that in most cases the binary system does not reach the Darwin instability during the 
AGB phase of the primary star, and so there is no CE phase at this stage. Namely, we consider 
systems where the secondary (initially less massive star) is massive enough to bring the primary 
envelope to synchronization with the orbital motion, and a CE phase before the secondary has 
developed a CO core is prevented . The primary forms a WD of mass M\ydi ^ O.65M . 

Later in the evolution the secondary evolves and interaction takes place only at late AGB 
phases. The WD remnant of the primary star cannot bring the secondary envelope to synchroniza- 
tion, and the binary system enters a CE phase. However, as this occurs late on the AGB when the 
core is a CO core and is massive enough, M core > O.6M . Despite the large uncertainties in the 
binary interaction, systems where the combined remnants mass is M\vdi +M corc > 1.4M are very 
likely to be descendant of two stars that each left the main sequence with a mass of Mms > 2.3M . 

On the other hand, WD binary systems with M\vDi+M core < 1.3M are likely to be descendant 
of systems where at least one of the two stars had a mass of Mms ^ 2.3M and a strong tidal 
interaction took place on the RGB. This would further prevent the WD from achieving a high mass. 



The WDs remnants, or a WD and the core of the secondary, of two stars in a binary system 
that are likely to have strong interaction on the AGB but not on the RGB, is about equal to the 
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Mprimordiol /Me 



Fig. 1. — The maximum radius Rr on the red giant branch (RGB) and the one Ra on asymptotic 
giant br anch (AGB) atta i ned b y starts as function of their main sequence (MS) mass, as calcu- 
lated by llben fc Tutukovl (|1985l ) . We added two thick dashed-dotted lines that are the quantity 
log(Rjy/Rii) as given by our approximate equations (pQ) and ([2]) for the two mass ranges, respec- 
tively. 
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Chandrasekhar mass limit. The condition is that both stars have Mms > 2.3M when they leave 
the MS, and the orbital separation be in the right range. The orbital separation range for forming 
a SN la progenitor is not simple to estimate, as it requires knowledge of the tidal interaction. In 
any case, in the next section we will use the logarithmic range A = \og{RjC) — log(i?i{) in estimating 
the number of systems that potentially can lead to SNe la. 

When a low mass secondary star, i.e., Mms < 2.3M , evolves, it has a lower envelope mass. 
A WD companion (the remnant of the primary) is more likely to survive the CE phase. The 
remnant is a close WD binary system with a total mass of Mj < 1.3M@. Those WD binaries with 
Mt > IAMq are likely to come from s ystems whe re the secondary had Mms > 2.3M , and the 
CE is more likely to end with a merger (|Sokerll2012l ). This is because the outcome of a CE phase is 
sensitive to the envelope mass. A small increase in the envelope mass above the combined ma sses of 
the c ore and WD companion can lead to a merger event at the termination of the CE phase (poker 
2012j). To have a core mass of M corc > 0.7 M the secondary MS mass should be Mms ^ 3M . 
The envelope mass is then ~ 2.3M which is 3 times and more the companion WD mass. 



This general discussion might explain the finding of iBadenes &: Maoa ([2012]) that the super 



Chandrasekhar merger rate of WD binary systems i n the s olar neighborhood , ~ 10 yr M 



is only ~ 10% of the SN la rate found bv iLi et all (|201lh . IBadenes fc Maozl (|2012h suggest that 



' 



the explanation is that sub-Chandrasekhar mergers lead to SNe la. To account for the observed 
rate they require SNe la to occur in merger of WD with a total mass of as low as ~ 1.1M . 
Ruiter et al.l (|201lh find in their population synthesis that theoretically there are potentially enough 
sub-Chandrasekhar progenitors to account for the rates of SNe la. Instead, we suggest that 
many single compact-magnetic super-Chandrasekhar WDs explode as SNe la. They are the rem- 
nants of core- WD merger at the termination of the CE phase. The recent population synthesis of 
Garcia-Berro et al.l (J2012J) show that many more mergers occur during the AGB phase than of two 
WDs (about a factor of 8). 



3. ASSUMPTION AND SET UP 



In this section we list the assumptions and set up of the population synthesis calculations. 
Because of the large uncertainties in the CE evolution and the controversial a-CE prescription 
(|Sokerll2012l ). we limit our goal to show that the number of SN la can be explai ned with the CD 



and DP scenario. Previous population synthesis of the the DD scenario (e.g., Yuneclson 



Ruiter et al 



(120121) 



and 



200 91; iMennekens et al.ll2010l ; IWang fc Hanll2012l ) are summarized bv iNelemans et al 



2005 



Wang &: Hani ([20120 . The more recent population synthesis study conducted by 



Toonen et al.l ([2012]) shows that the expected total number SN la progenitors in the DD scenario 
is a factor of ~ 7 — 12 times lower than observed. Even if the fraction of binary is taken to be 
70% instead of 50%, th e contribution of t he DP scenario is < 20%. This is compatible with the 
observational finding of IBadenes &; Maoa (|2012l ). With the CE prescription used in the context 
of the PP scenario, and for 50% of systems in binaries, the results span the range of 0.2 — 0.6 SN 
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Ia per IOOOMq star formation, which is below half the observed value. As we will show below, the 
CD scenario does much better and can account for the observed value. 



There a re three basic differe nces between our calculations and those summarized bv lNelemans et al 



(|2012l ) and IWang &: Hani (|2012l ). (i) We consider the merger of the WD companion with the core, 
i n particular in cases where the WD mass is less than the core mass, to be of very high probability 
(|Sokerll2012l ). (ii) We use a simple prescription for the relevant binary orbital separation that might 
allow stronger interaction on the AGB (see more in section [5]). (in) Our usage of a high efficiency 
mass transfer (77 > 0.5 below) is more than what other population synthesis codes use/find. This 
effect might be the lar gest of the thre e effec ts listed here. In that regards we note that our results 
disagree with those of iMeng &: Yangl (|2012l ) regarding the CD scenario, both quantitatively and 
qualitatively. It is possible that most of the differences result from our usage of this efficient mass 
transfer process. 



3.1. Binary evolution 



The ingredients of the population synthesis co de are as follow s. 



(1) Initial mass functio (IMF). We use the IMF of iKroupa et al.1 (|1993l ) where 



f(m) 



0.035m- L3 , for 0.08 < m < O.5M 
0.019m- 2 - 2 , for 0.5 < m < l.OM 
0.019m- 2 - 7 , for 1.0 < m < 00, 



(3) 



where m = M/M^ 
Kroupa et all (|l993h 



In our numerical monte-carlo code we take the formula as suggested by 
(0.19X 1 - 55 + 0.050X - 6 ) 



M = °- 08 + (1 _ X) 0.58 

where X is a uniformly distributed random number in the interval [0, 1]. 



(4) 



(2) Upper mass limit of M < 7M & . From these stars we choose the ones in the mass range 
2.3Mq < M\ < 7Mq. These are the primary stars in what follows. The secondary in each binary 
system is taken to have an equal probability to be in the mass range 0.08M Q < Mi < M\. These 
masses take into consideration that primaries with a mass of M\ > 7M & either explode as a core 
collapse supernova (CCSN) or form a NeMgO WD. 

(3) Lower mass limits. Once M\ evolves it might interact with the companion. To form a CO 
WD we require that a strong tidal interaction takes place on the AGB phase of the both stars, 
but not during the RGB phase. For the statistical purposes of this study we are not interested in 
the value of the orbital separation where tidal interaction takes place, but only in the ratio of that 
orbital interaction on the AGB to that on the RGB (see section [2]). Crudely, this ratio is the ratio 
of the maximu m radius of the s t ar on the AGB, Ra, to that on the RGB, Rr, as given by equation 
([2]) taken from llben fc Tutukovl (|1985l ). For M < 2.35M Q this ratio is very close to 1, and we ignore 



- 8 - 



systems with Mi < 2.3M . After mass transfer takes place we make sure that M2 new > 2.3M . 
We ignore systems with low initial secondary mass of M 2 < 1M . Low mass secondary stars that 
are potentially in the relevant orbital separation will not be able to bring the primary envelope to 
synchronization with the orbital motion, and will enter a CE phase too early. 

(4) Mass transfer. The strong interaction during the AGB will lead to mass transfer from the 
primary to the secondary. We consider here a channel where the mass transfer from the primary 
to the secondary avoids a CE phase. A CE phase occurs only when the secondary becomes a giant. 
Mennekens et al.l (|2010l ) find this channel to be more common than the two-CE phases channel in 
their population synthesis study of the DD scenario (some details and examples of the evolutionary 
channels can be found there). We parameterize the poorly known mass transfer process as follows. 
We take the new mass of the secondary after the primary has gone through the AGB phase and 
turned into a WD of mass Mwdi as 



M 2new = M 2 + r j (M 1 - Mwdi] 



(5) 



where < r? < 1 . Fro m a more accurate population synthesis, using the code described by 
Garcia-Berro et al.l (|2012l ). of the mass transfer process from the evolved primary to the secondary, 
the average value of rj is found to be fj ~ 0.9 (Camacho, J. 2012, private communication). We will 
take rj = 0.8 as our standard case. 

(5) WD Masses. The mass o f the WD that emerges from the evolution of each star is taken 
according to ICatalan et al.l (J2008J) as 



Mwdi 



(0.096 ± 0.005)Mi + 0.429 ± 0.015, M 1 < 2.7M 
(0.137 ± 0.007)Mi + 0.318 ± 0.018, M x > 2.7M , 



where for the WD remnant of the secondary M\ is replaced by M 



2new • 



(6) Binary fraction. Our standard fraction of binaries for these massive stars is taken to be in 
the range fb = . 5 — .65, where this fraction includes secondaries in the mass range O.O8M — M\. 
Nelemans et al.l (J2012J) use fs = 0.5, but our condition of M\ > 2.3M and M 2 > 1M takes more 
massive stars, where the binary fracti on can be somewhat la rger that = 0.5, up to ~ 70% (see 



summary of observations in fig. 12 of lRaghavan et al.1 (|2010l ). 



3.2. Counting SNIa Progenitors 



The criteria to count a system as a potential SN la progenitor are as follows. 
(1) Mass considerations. For the SN la progenitors we count only systems with 

M WD i +M core > 1.4M . 



(6) 



In most works of this type th e final orbit al separation is calculated using the controversial CE 
a-prescription. We here follow ISokerl (|2012l ) and assume that most cases where the AGB envelope 
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is massive enough will end with merger. This is particularly the case when the WD remnant of 
the primary star has a lower mass than the core of the secondary, namely, Mwdi < M corc . These 
systems will be counted separately as the likely progenitors of SN la in the CD scenario. 

(2) Orbital separation considerations. We consider potential progenitors of WD-core or WD- 
WD merger to be systems (namely, systems where the WD enters the envelope of the secondary 
during the AGB phase of the secondary star) that have strong tidal interaction on the AGB but 
not on the RGB. The relevant initial orbital separation will be in the range C q Rr < a < C q Ra, 
where Rr and Ra are the maximum radius of the star on the RGB and AGB, respectively. The 
factor C ~ 2 — 5 is determined by the tidal forces. The relevant range for our potential progenitors 
in the logarithm of the orbital separation scale is therefore log(C i?A) — log(C i?jj) = \og{RA/ Rr)- 
We take the initial orbital separation of the binary population of massive stars created by the 
numerical code to span a range of 4.5 orders of magnitude (from a m \ n ~ 10i? Q to a max ~ 1500 AU) 
with an equal probability in the logarithmic of the orbital separation. Accordingly, the probability 
of a binary system to have the initial relevant orbital separation is 



Xi 



\og{C R A ) - log(C R R ) 1 R A I 

~ — log = (2.2 

log a m ax - log a m in 4.5 R R 4.5 



1.8 log AW) 



(7) 



where equation [2] based on llben fe Tutukovi (jl985l ) have been used in the last equality. Instead of 
calculating for each of the two stars, in the present study it is adequate to take the average mass 
for this criterion M sta r = 0.5(Mi + M 2 ). 

(3) Progenitors of CCSNe. Systems that have M\ > 8.5M^ are co unted as having one CCSN, 
while if M2 > 8.5M0, the system produces two CCSNe (jGilmore 1120041 ). This number will be used 
to find the ratio between the nu mber of SN la and that of CCSN, and compare with observations, 
e.g., as given bv lLi et al.l (|201ll ). 



(4) Avoiding NeMgO WDs. Merger remnants that have a large fraction of neon and magnesium 
cannot explode. We therefore avoid NeMgO WDs by not counting systems where either Mwdi > 
1.1M or M C ore > 1-IMq. 



4. RESULTS 

4.1. The number of potential SN la progenitors 

We consider the parameters range of the model to be /& — 0.5 — 0.65 and r] ~ 0.5 — 0.9, for 
the fraction of binary systems and mass transfer parameter, respectively. The number of systems 
obtained under the counting conditions listed in section 13.21 are presented in Table [TJ The first 
panel gives the number of potential progenitors of SN la per 1000M Q of stars formed, while the 
second panel gives the number of potential progenitors relative to the number of CCSNe. We also 
give the numbers only for systems that in addition to the criteria given in section [31 also obey the 
condition that the WD mass is lower than the core mass, Mwdi < Af core . Such systems are more 
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likely to merge at the termination of the CE phase. 

The number of potential progenitors for the standard parameters and with the condition of 
-Mwdi < -^corei is ~ 1 — 2 times the observed value. There are several processes that can change 
the number of potential SN la progenitors in our modelling, but not by much, (i) Some cases 
with Mwm > M cme can also lead to SN la, although not many such progenitors are expected 



(poker et al.ll2012l ). On the other hand, we might overestimate the number of progenitors as during 
the merger process some mass might be lost, hence some systems will become sub-critical in their 
final merger product mass. As well, in some cases the primary WD enters the envelope before the 
core reaches its maximum mass on the AGB, M core . Hence, some of the systems counted here as 
having a total mass above the Chandrasekhar mass, might actually be below the critical mass, (ii) 
The binary fraction can be larger than ft, = 0.5, an d /r ~ 0.65 migh t be m ore appropriate for the 



B-type stars considered here (e.g., see summary bv lRaghavan et al.1 (|2010l )). On the other hand, 
in some systems a tertiary star might disturb the evolution and prevents a SN la from occurring. 
(to) A value of r/ = 0.9 is more appropriate. This will increase the number of progenitors. On 
the other hand, some unstable phases in binary evolution might lead to extra mass loss that is not 
considered in population synthesis calculations. As well, some systems might enter a CE evolution 
during the AGB phase of the primary. This might prevents the secondary star from developing a 
CO core later on. 

Overall, we consider it very satisfactory to the CD scenario that with a very simple modelling 
we find the number of SN la to be slightly above the observed value. This margin comes for the 
over-counting mentioned above. More sophisticated population synthesis calculations that take 
mergers into account should be performed to derive more accurate values for the expected number 
of SN la in the CD scenario. As well, the time delay should be calculated. Here we only compare 
the total number (time integrated) of SN la, and not the evolution of rate with time. The model 
should eventually reproduce the delay ti me distribution of SN la if it to account for most SN la. 



If the merging core- WD events found by iMennekens et al.l (|2010l ) did not consider the many core 



WD merger events they found to be SN la. If these events are considered to be SN la, then their 
calculations yield a similar rate to the one we obtain here under the same conditions (N. Mennekens 
2012, private communication). 



4.2. Properties of the merging core- WD systems 

We present the distribution in some properties of the WD-core remnants. All graphs in this 
subsection are for systems that obey the conditions of section 13.21 and for the mass transfer param- 
eter of r\ = 0.8. In these graphs 500 systems correspond to 0.073 systems per 1000M Q of star 
formation for /b = 1 (all stars are binary). For other values of /b the number should be scaled 
accordingly. 

In Figure [2] we present the mass histogram of the total mass Mt = Mwdi + -^core for systems 



Table 1: The calculated number of systems entering the CE phase according to the counting criteria 
listed in secti on 13.21 Upper pane l: number of systems per IOOOMq of star formation. The observed 
value is from iMaoz et al.l (|2012l ). Second panel: The calculated ratio of the number of p o tentia l 
SN la p r ogeni tors to core collapse SNe (CCSN). The observed range is from IMaoz et al.l (|201ll ). 
Li et al.l (|201ll ) and iMelinder et al.l (J2012J). The two rows labelled with superscript * in both panels 
are for systems with il%Di < M COTe , that are more likely to suffer core- WD merger and have long 
delay to explosion. 



SNIa per 1OOOM (Observed value 1 - 2) 
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rj = 0.8 
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SNIa/CCSN (Observed value 0.2 - 0.4) 
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T) = 0.8 
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0.44 


0.51 
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0.26 


0.44 


fb = 0.65 


0.55 


0.65 


fb = 0.65* 


0.32 


0.55 
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with Mwdi < M core . In most of these systems the cold primary WD is expected to be destructed 
on the core. A long ti me delay to explosion i s likely. We emph asize that the merger process might 
eject some mass (e.g., lHan fc Webbinklll999l ; IChen et al,ll2012j ). and the final merger product mass 
might be smaller. Due to mass loss during the merger process, we expect that more systems will 
end in the mass range 1.4 — 1.5M@, where they expect to have a very long time delay to explosion 
(up to ~ 10 10 yr), long after the star formation episode ceased. Systems with mass much above 
the Chandrasekhar mass limit, but where the WD is destructed on the core, will have shorter 
time delay, but longer that the time the ejected envelope disperses in the ISM, ~ 10 5 yr, and are 
expected to explode in star-forming regions. Systems where Mwdi ^ M core are less likely to merge. 
If they do merge, then the core might be destructed on the cool er WD, and explo sion might occur 
within tens of years, leading to the presence of hydrogen lines (jSoker et al.ll2012l ). The histogram 
of the mass ratio for any remnants mass ratio is given in Figure EJ 

In FigureHJwe show the density of potential SN la progenitors in our model in the Mwdi — M corc 
plane. In Figure[5]the density in the initial mass of the binary components plane is given for systems 
where Mwdi < M core . By definition, all systems have Mi > Mi- However, in our modelling of 
the CD scenario most, but definitely not all, potential SN la progenitors have Mwdi < M core . As 
expected, the core mass is not much larger than the primary WD mass in the merging systems. 

In Figure [6] we show the ratio of the envelope mass to the WD mass, £ = (M2 new — M core )/MwDij 
as function of core mass, for cases with Mwdi < M core (upper panel) and Mwdi > M core (lower 
panel). Neglecting mass loss by the secondary star before the onset of the CE phase, £ is the ratio 
of secondary envelope mass t o the prima ry stellar remnant mass. High values of £ > 3 — 4 favor 
merger during the CE phase (|Sokerll2012l ). It is evident that systems for which Mwdi < M core are 



more likely to merge. That the hot core is more massive than the cold WD companion makes a 
difference in the merger process relative to the process of merger of two cold WDs. 



5. DISCUSSION AND SUMMARY 

In estimating the number of SN la progenitors in the core-degenerate (CD) scenario we took 
a different approach than previous population synthesis studies in treating the common envelope 
(CE) phase. We assumed that most CE processes where the mass of the white dwarf remnant of 
the primary (initially more massive) star is less that the mass of the core, Mwdi < M core , end in 
mergers. These systems have a large envelope to WD mass ratio, £ > 3.5, as can b e seen in th e 



upper panel of Figure [6) Such systems have a high probability for WD-core merger ((SokerJ |2012j) . 
This approach avoids using the highly uncertain and controversial a-CE parameter in estimating 
the final orbital separation. 

We required that the binary interaction is strong during the late AGB phase of the secondary 
star, but not during its RGB phase. We then used the ratio of maximum radii on these two 
phases, as presented in Figure [IJ to find the probability for each binary system to be in the 
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relevant orbital separation range. The probability is given in equation (|7|), where we took the 
initial orbital separation of the binary system population to span 4.5 orders of magnitudes. This 
simple treatment avoids using the tidal interaction expression that has some uncertainties. For 
having a large ratio of AGB maximum radius to RGB maximum radius the initial mass of the 
star should be Mms > 2.3Mq. Such stars also form WDs of masses that are about half the 
Chandrasekhar critical mass, such that the merger product of two such stars reaches the critical 
mass. This coincidence is discussed in section [2j 

These steps, that are described in section [31 leave us with two parameters. One is the mass 
transfer parameter from the primary to the secondary, 77, as given by equation [SJ and the other is 
the binary fraction The likely values for these are r/ ~ 0.9 and ft, ~ 0.65. We ta ke a conservative 



appro ach and use 77 = 0.8 and fb = 0.5 — 0.65. The value of fs = 0.5 is used by iNelemans et al 



([2012) in their comparison of different studies. The criteria for a system to be counted as a SN la 



progenitor are listed in section 13.21 

The results are given in Table [TJ The rows without a superscript * include all core- WD 
systems entering a CE phase, while the rows marked with a superscript * include only systems 
with Mwdi < M core . The latter cases are more likely to merge and are the systems considered 
here. These numbers should be taken very cautiously at this preliminary stage. As can be seen from 
Table [H our very simple approach matches rate of SN la deduced from observations. In section I4TT1 
we also list processes we have neglected, some that will increase and other that will decrease the 
numbers found here. Overall, despite the large uncertainties and crude estimate, we find it very 
satisfactory that the simplest population synthesis we could think of yields SN la average rate in the 
CD scenario very close to observations. This match must be reexamined with more sophisticated 
population synthesis calculations that take mergers into account and where onset of the CE phase 
is reconsidered. 



Nelemans et al.1 (|2012i ) summarize the results of population synthesis studies of the double- 
degenerate (DD) scenario from six groups. The results span the range of 0.2 — 0.6 SN la per 
lOOOM© for fs = 0.5. These numbers are a factor of 2-10 smaller than our results. We think that 
two processes combine to give these lower numbers. First is the CE-a prescription they use. Some 
systems in their calculations end at a too large orbital separation to become a SN la. The second 
one has to do with the treatment of the tidal interaction and mass transfer during the RGB and 
AGB phases. The tidal interaction and the formation of a CE phase when the secondary reach 
the upper AGB might be more significant than what is usually assumed in population synthesis 
studies. This is because of the unstable nature of massive AGB stars, as all secondaries are in our 
calculations. Such instabilities are likely to increase the number of WD that enter the CE phase 
with the secondary. We call for people who follow the RGB and AGB phases in the simulations 
to reconsider the tidal interaction when the AGB reaches the upper AGB, and alow for stronger 
interaction due to instabilities. Such instabilities might increase the orbital separation from where 
the WD remnant of the primary star enters the CE phase. We intend to study the effects of AGB 
instabilities in a future paper. 
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Future studies of the CD scenario will have to explore the exact nature and outcome of the 
merger process of a hot core with a cooler WD with a similar mass 0.7M core < Mwdi ^ l-2M core , 
as this is a major open question in the CD scenario. Three-dimensional simulations of this process 
are highly desired before the CD scenario can stand on a more solid ground. 

On the observational side, if the CD scenario does indeed account for a substantial fraction 
of SN la, then rapidly rotating super-Chandrasekhar WDs should be found. Claims for massive 
rapidly rotatin g and stro ngly magnetized WDs have already been made, e.g. see discussion by 
Malheiro et al.l (|2012l ) and IWickramasinghe Ferrarid (120001). Out of the 16 magnetic WDs with 
well determined mass reported by IWickramasinghe &: Ferrarid (|2000i ). two have masses close to 
the Chandrasekhar limit, M > 1.3Mq. Recent st udies also show that magnet ic WDs on average 
are more massiv e than non-magnetic WD s (e.g., IVanlandingham et al.1 120051 ). in particular for 
-M\VD > 1-3M Q (jNalezvtv fc Madejl 12004 ). WD with a long time delay to explosion should be 
found in the mass range ~ 1.4 — 1.5M Q (Jorge Rueda, 2012, private communication). More 
massive merger remnants are expected to explod e in star forming r egions, and in some cases even 
have massive hydrogen-rich nebula aro und them (poker et al.ll2012i ). For further discussion on the 
expected fraction of massive WDs see lllkov Sokeri (|2012h . There are also theoretical studies 
of su per-Chandrasekhar WDs stabilized by strong magnetic field (e.g., iKundu Mukhopadhvav 
2012J) that should be extended to systems relevant to the CD scenario. 
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1600 




^WD.+^core 

Fig. 2. — Histogram of the combined masses of the two remnants M± = Mwdi + M core . Only 
systems with Mwdi < M core are included in this figure. On average, these are more likely to merge 
at the end of the CE phase. 
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Fig. 3. — Histogram of the WD to core mass ratio, MwDiA^ccnre for the systems studied in this 
section. Note that a large number of systems have the AGB core mass larger than the mass of the 
WD during the final CE phase. 
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Fig. 4. — Density of number of systems in the Mwdi — -^core plane. The scale is the number of 
systems in each square (of 0.025 x 0.025Mq). Only in this figure we have systems with Mwdi > 
I.XMq and/or M core > 1.1M® included. This is for comparison purposes only, as systems where 
either Mwdi > I.IMq or M core > 1.1M are not counted as SN la progenitors. 



- 21 - 




Fig. 5. — Density of number of systems in the initial masses plane (the Mi — M2 plane). The scale 
is the number of systems in each square. 
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Fig. 6. — The ratio of secondary envelope mass to primary WD mass, ^ = (M2 new — M core )/M\vDi, 
as function of M core for rj = 0.8. Upper panel: As with most figures in this subsection, only systems 
with Mwdi < M core are included. Lower panel: Only syst ems with M wm > M core are included. 
High values of £ > 3 — 4 favor merger during the CE phase (jSokerll2012n , M2 new is the mass of the 
secondary after the mass transfer episode from the primary star according to equation ([5]). 



